Spontaneous trimerization in two-dimensional antiferromagnets 



o 
o 

(N 
D 

\o 



I 

o 
o 



> 
m 
(N 
cn 

ON 

o 

OO 

o 



Zi Cai, Shu Chen, and Yupeng Wang 
Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, 
Chinese Academy of Sciences, Beijing 100080, P. R. China 
(Dated: Received February 26, 2009) 

In this paper, we propose a novel quantum paramagnetic state in two-dimensional antiferromag- 
nets: the spontaneously trimer state, which is constructed by the trimers of spins. Each trimer is a 
singlet state formed by three neighboring spins with SU(3) symmetry. A frustrated spin-1 Heisen- 
berg model in the kagome lattice is investigated. By analogy to the pseudo-potential approach in the 
fractional quantum Hall effect (FQHE), we find that the trimer state provides a perfect description 
for the exact ground state of this model. Other interesting properties, such as the local excitations 
as well as magnetization plateaus have also been investigated. 

PACS numbers: 75.10.Jm, 71.27.+a, 75.10.-b 



The nature of quantum paramagnetic phases of 
the two-dimensional (2D) antiferromagnetic systems at- 
tracted considerable attention in the past twenty years 
because of its potential association with the physics of 
the cuprate superconductors [H. Because of absence of 
long-range order in the quantum paramagnetic phase, a 
central problem is the classification of the phases and 
critical points. According to their symmetries, we could 
classify these phases into two classes. The first kind is 
known as "spin liquid state", which restores the SU(2) 
symmetry of the Neel state and does not break any sym- 
metry of the original Hamiltonian. One of the well-known 
state in this class is the resonating valence bond (RVB) 
state or the "spin liquid" state [1]. 

On the other hand. Read and Sachdev have presented 
another possibility 1?| based on the Schwinger boson anal- 
ysis on the SU(N) quantum antiferromagnets. Their key 
point is that the condensation of the instantons with the 
Berry phase leads to the spontaneous breaking of the 
spatial symmetry (translational or rotational symmetry) 
of the original Hamiltonian. Similar with the solid, this 
state possesses a short range order and thus is nominated 
as valence-bond solids (VBS) state. As shown in Ref.[4], 
the structure of the VBS state depends on 2S(mod 4), 
where S is the spin of the SU(N) model. For the usual 
spin 1/2 systems, corresponding VBS state is known as 
"dimer state" , where two nearest spins form a singlet or a 
dimer: {][ — |t)/\/2 and the overall state is the product 
of all those dimers. 

The dimer state is first proposed as the exact ground 
state for the Majumdar-Ghosh(MG) model ^Si], which is 
an antiferromagnetic spin chain with nearest and next- 
nearest neighbor interactions. For a particular ratio of 
those exchange interactions, the model has a two-foldly 
degenerate dimer ground state. Inspired by the exact 
solution of the MG model, Shastry and Sutherland pro- 
posed a two dimensional model known as Shastry- 
Sutherland (SS) model, whose ground state is the exact 
dimer state in 2D. It has been used in understanding the 
physical properties of SrCu^{B0^)2 01j which is topo- 
logically equivalent to the SS model. 

Most above researches focus on the spin-1/2 system. 



because of its potential correlation with the high-Tc su- 
perconductivity. However, the cold atom systems have 
provided an ideal playground to investigate the high-spin 
strongly correlated system experimentally. For example, 
several proposals have been provided to realize the spin-1 
Heisenberg model in the optical lattice [8j]. Even in one 
dimension, new possibility emerges and a typical exam- 
ple is known as the Affleck-Kennedy-Lieb-Tasaki (AKLT) 
state 9] . This state breaks no symmetry of the Hamilto- 
nian and is believed to provide a good description for the 
ground state of the pure spin-1 Heisenberg chain [13, EH- 
As to the 2D case, the situation is apparently more com- 
plex, the VBS state in the spin-1 quantum paramagnetic 
phase could be either the dimer state [l2| that sponta- 
neously breaks the translational symmetry, or a gener- 
alized AKLT state [ll| that only breaks the rotational 
symmetry. 

In this paper, we present a novel VBS state named as 
"trimer state" for the spin-1 system in the 2D Kagome 
lattice. Such a state is a product of trimer singlets with 
each spin singlet formed by three spins with the total spin 
zero. Inspiring by the seminal work of Ref. [l^, we con- 
structed a frustrated Heisenberg model in a kagome lat- 
tice, and found that at a particular ratio of the different 
couplings, the "trimer state" provides a perfect approxi- 
mation for the exact ground state of this model. Similar 
states have been discussed in an SU(3) spin tetrahedron 
chain [3, H, [l 5| . [l ^ and in the kagome lattice with dis- 
torted coupling[l7|. To our knowledge, we provide the 
first example in an isotropic 2D SU(2) antiferromagnets 
with uniform coupling, where the spins are spontaneously 
trimerized. Then we investigated the excited properties 
and the magnetization plateaus in the applied magnetic 
field. It is believed that this state represents a novel class 
of the VBS state: N neighboring spins cluster together to 
form a singlet with SU(N) symmetry and spontaneously 
break the spacial symmetry. Within a cluster, each spin 
is maximally entangled with other spins. We call this 
kind of state as "cluster state" while the usual dimer 
state is the simplest case. 

The ground state of the quantum antiferromagnetic 
Heisenberg model on the kagome lattice has been in- 
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FIG. 1: (a) The structure of frustrated Heisenberg model 
in kagome lattice, with the nearest coupling J, the second 
nearest coupling Ji and the third nearest coupling J2. (b) 
One block consists of 5 sites and each one equally couples 
with others 

tensively investigated by various methods, but is still 
far from being totally understood even for the spin-1/2 
case. Recently, Laws et. al found a A^i^+-based material: 
Ni3V20s [1^, which is topologically equivalent to a frus- 
trated spin-1 Heisenberg model in kagome lattice. The 
spin-1 model on the kagome lattice has been investigated 
in several works [l9l I20I [2ll . [23 | , but most of them focus 
on the uniaxial anisotropic case. Here we focus on the 
SU(2) case and introduce the second nearest ((ij)) and 
the third nearest [ij] couplings Ji and J2. Explicitly, the 
Hamiltonian for this model is given by 

i/ = J^Si-Sj + Ji ^ -S. + J2^Si-Sj, (1) 

(iJ> ((iJ» [iJ] 

with the structure schematically shown in Fig.l. 

In this paper we would focus on the point J = 3 Ji — 
3 J2 and investigate the ground state at this point by the 
"pseudopotential approach" , which was first used to get 
the famous Laughlin wavefunction in the FQHE [2^ . The 
analogy between the Heisenberg model and the FQHE 
was first introduced by Arovas et. al. in the ID 

case, where the Hamiltonian of a spin-1 antiferromagnetic 
Heisenberg model is decomposed into the summation of 
projection operators 

=^Si.Si+i -^[3p2(j,i + i) + pi(j,i + i)_2]. 

i i 

(2) 

Observing that the AKLT state [9] is the exact ground 
state of the first part of the summation of projection 
operators J^i ^P'^ih * + l)i they take it as a trial ground 
state of the spin-1 Heisenberg chain and consider the sec- 
ond part J2i P^ih i -f 1) as a perturbation. It turns out 
that the AKLT state is a very good approximation of 
the exact ground state of Eq. (2) and numerical results 
show that the difference between the ground state en- 
ergies of these two states is within 5%. Recently, we 



generalized this method to the 2D spin-1 Ji — J2 an- 
tiferromagnetic model and found that at the maxi- 
mal frustrated point (Ji = 2J2), the ground state could 
be described by a two- fold 2D generalized AKLT states, 
which completely agrees with the general prediction by 
the field theo ry[4| , and was verified by the numerical re- 
sults recently |24|. 

Next we use the similar method to investigate the 
ground state of Eq. (1). Notice that at the point 
J = 3Ji = 3J2, Eq. (1) could be rewritten as the sum of 
identical blocks 58q, as shown in Fig l.(b). Each block is 
constructed by 5 spins and every spin is coupled to the 
other four spins identically, thus we have 

where a denotes the ath block. Then we expand this 
block Hamiltonian by the projection operators of the to- 
tal spin in a block. In terms of the projection operators, 
we have 

5 

*B„ = X CsPl - 5 with Cs = S{S + l)/2. (4) 

s=o 

To make progress, we further decompose it as two part 

with »o = 15P5+10P^ + 6P^-5and »i =3P2+pi. 
The operator projects the spin state of the ath Block 
onto the subspace with total spin S. Dividing the origi- 
nal Hamiltonian into two parts is inspired by the success 
of Haldane's pseudopotential method in FQHE ^2^ and 
the model of the spin-1 chain 10]. Return to our spin 
model, now it is clear why we divided the Hamiltonian 
like this: notice that the coefficient decreases rapidly as 
S descends, thus if we can find the exact ground state of 
So = J2a^a, we can treat the left part »i = J2a^i 
and investigate the properties of the ground state. 
First we focus on the exact ground state of the model 

»o = E[l5P^ + 10P^ + 6Pa-5]- (5) 

a 

Since *8|^ is positive semidefinite, any state with the total 
spin of each block S'J < 2 is the exact ground state of 
Eq.(5). It is not difficult to find that the only possibilities 
are the two-fold degenerate states as shown in Fig. 2(a) 
and (b), where each block possesses a trimer singlet. Ex- 
plicitly, the trimer singlet can be represented as 

which is a singlet composed of three spins on neighboring 
sites i, j, and k with total spin zero (Fig. 2(c). Here 
denotes the spin on site i with the value a taking 1, 0, or 
— 1 and ea[3-y is an antisymmetric tensor. Because there 
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FIG. 2: (color online) (a) and (b): The two-fold degenerate 
ground states: \'i'a) and I'I't). (c) The structure of one trimer: 
a AKLT chain with Length L=3 and periodical boundary con- 
dition 

exists a trimer singlet in each block, the total spin of 
each spin block could not be larger than 2. Therefore it 
is straightforward that we have found the exact ground 
states of *8o, or in other words, the trial ground state of 
our original Hamiltonian Eq. (3). 

Now we will discuss the effect of the perturbation part. 
Since the trial ground states are two-fold degenerate, 
there seems to be a possibility that the perturbation 
would resonate these two trial ground states jvPa) and 
Iv^b) to form a linear superposition and further lower the 
energy of the ground state. However, we would show 
this is not the case, at least in the thermodynamic limit. 
Let's calculate the iroirdiagonal term {'i>a\H\'^i,), it is not 
difficult to find that this nondiagonal term vanishes in 
the thermodynamic limit(A^ oo), just like the ID MG 
model. So the perturbation would not shift the two-fold 
degeneracy of the trial ground states and we can safely 
conclude that \'^a) and provide a good description 
for the Hamiltonian Eq. (3). It is straightforward to get 
the ground energy of this variational ground state at the 
point J = 3 Ji = 3 J2 

{^a\H\^a)=-ljN, 

where N is the number of the triangles in the lattice. 

It is also interesting to study the properties of the 
excited states as well as magnetization process of this 
model. It's well known that for the 2D dimer model, 
there are magiretization plateaus in the magnetization 
curve because of the localization of the single triplet 
excitation 0, H^. It is natural to ask whether it 
happens in our model. Notice that for a single trimer, 
the lowest excitation state is the triplet state, which is 
the product of a dimer singlet composed of two spins 
(-i=(|l,-l) + I - 1,1) - |00)) and a single spin state 

(|1),|-1) or|0)). 




FIG. 3: The triplet excited states: a triplet locates at m (a) 
and n (b), m and n are nearest neighbor triplets. A triplet is 
formed by a dimer of two spin 1 and a single spinon (c) 

The triplet state is shown in Fig. 3. For convenience, 
we use I to) denoting the TOth trimer singlet being excited 
to a triplet. To study the low energy excitation of this 
model. We calculate and find that 

{m\H\n) = (6) 

for m ^ n and 

{m\H\m) = -J-p{N-l) (7) 

where H is the Hamiltonian (3), and |to) and \n) are 
shown in Fig. 3 (a), (b). Eq.(6) means that the single 
triplet excitation is almost dispersionless, at least to the 
order we considered. Eq.(7) means that it is gapped, 
with a gap of ^ J. This dispersionless triplet excitation 
in our model is different from its one dimensional ana- 
logue [13, [11], and it would localize or form a bound 
pai r with other single triplet, just as in the dimer state 
poL [30I I . Because of the localization of the single triplet 
excitation, there are magnetization plateaus in the mag- 
netization curve in our model. The naive picture for the 
magnetization process is shown in Fig. 4. Therefore, at 
least two magnetization plateaus appear at TO/TOsat=l/2 
and 3/4 corresponding to the phase (b) and (c) in Fig. 4. 
It is known that similar magnetization process and mag- 
netization plateaus have been studied in other Spin-1 
systems, such as S=l spin chain [28j. or S=l Heisenberg 
model (only including the nearest coupling) in uniform 
and distorted kagome lattice 31]. 

However, we have not considered the effect of the in- 
teractions between the triplets, which may lead to new 
plateaus. Further more, it is possible that the higher or- 
der interaction would result in the pair of neighboring 
triplets, similar to the dimer phase (29l. [soj. 

Now we will discuss the stability of our trimer state. 
Up to now, all of our analysis are based on a special point 
of our coupling parameters: . At this point, the trimer 
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FIG. 4: The magnetization process for a single trimer.(a)The 
ground state with total spin Stotai=0 .(b)The state with to- 
tal spin Stotai~i, it is formed by a dimer and a spinon with 
S=l(the bigger red arrowhead). {c)Stotai='2, the smaller blue 
arrowhead represents the spinon with S=l/2. {d)Stotai=3, 
the state with saturated magnetization 



biquadratic term, actually provides a perfect variational 
ground state for the pure S=l Heisenberg Hamiltonian 
without biquadratic term. Therefore in our case, the per- 
turbation of including the small spin biquadratic term 
would not change the nature of our trimer state. 

la) -w . lb) 




state actually provide a good approximation of the exact 
ground state, and slight perturbation should not change 
the nature of this ground state because it is a gapped 
VBS state which can not change by small perturbation 
unless it can overcome the energy gap, just as in one 
dimensional AKLT state, introducing small next nearest 
coupling can't change the nature of ID VBS statefll] . 

Usually there are two kinds of perturbation, the first 
kind is deviation from the special point of the cou- 
pling parameters J — 3Ji — 3J2, as analyzed above, 
where slight deviation will not change the nature of the 
ground state. However, there is another important point: 
Ji = J2 = 0, which corresponds to the S=l antiferromag- 
netic Heisenberg model with only nearest coupling. Hida 
first studied the ground state in this point by means of 
exact diagonalization and the cluster expansion [T^. It 
is shown that the ground state in this case is the hexagon 
singlet solid (HSS) state, rather than the trimer state in 
our case. The second possible perturbation for this kind 
of spin-1 system is the spin biquadratic term, which is 
necessary to construct the ID AKLT modelQ. However, 
as shown in Ref . [T^ , the exact ground state of the AKLT 
Hamiltonian {3P^(i,i + 1) in Eq.(2)) which includes the 



FIG. 5: (a) The cluster state (N=4) in spin-3/2 ladder. 
(b)The structure of an SU(4) cluster singlet. A 3/2 spin can 
be considered as a symmetrized state of three 1/2 spins. The 
bold lines represent a singlet for two 1/2 spins 



The trimer state is not the only example of our cluster 
state, another example is the spin-3/2 SU(4) spin ladder 
[33! . At particular value of the coupling, the ground state 
of this model is the exact spontaneous plaquette ground 
states, in which four 3/2 spins form an SU(4) spin singlet 
plaquette and spontaneously breaking the translational 
symmetry, as shown in Fig. 5. This model is a generaliza- 
tion of the MG model to the spin ladder systems and its 
ground state belongs to our cluster state with iV = 4. 

After completing this work, we learnt that in a recent 
paper 3J| , similar topic about the 2D multipartite VBS 
state has been discussed. Z.C. thanks Flavio S. Nogueira 
for helpful discussions. This work is partially supported 
by NSF of China under Grant No. 10574150, MOST 
grant 2006CB921300 and programs of Chinese Academy 
of Sciences and Deutscher Akademischer Austausch Di- 
enst (DA AD). 
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